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We show that the pinning of collective charge and spin modes by impurities in the cuprate 
superconductors leads to qualitatively different fingerprints in the local density of states (LDOS). In 
particular, in a pinned (static) spin droplet, the creation of a resonant impurity state is suppressed, 
the spin-resolved LDOS exhibits a characteristic spatial pattern, and the LDOS undergoes significant 
changes with increasing magnetic field. Since all of these fingerprints are absent in a charge droplet, 
impurities are a new probe for identifying the nature and relative strength of collective modes. 

PACS numbers: 74.72.-h,74.25.Jb,73.20.Mf, 74.25. Ha 



Whether the unconventional properties of the high- 
temperature superconductors (HTSC) arise from the in- 
teraction of electronic degrees of freedom with collective 
spin modes, charge modes, or phonons, is one of the 
key issues in understanding these complex materials 
While some angle-resolved photoemission (ARPES) ex- 
periments observed effects of the magnetic resonance 
mode on electronic excitations, other ARPES stud- 
ies argued that the electronic dispersion exhibits a num- 
ber of prominent features that arise from the coupling to 
phonons 0. Moreover, recent scanning tunneling spec- 
troscopy (STS) experiments 6] have reported evidence 
for a collective mode whose frequency is doping indepen- 
dent and shifts with isotope substitution, two properties 
that are inconsistent with those of the magnetic reso- 
nance mode 0, 0|. The interpretation of these exper- 
iments is further complicated by the coupling between 
charge modes and phonons [3, S|- Clearly, new experi- 
ments are required that can unambiguously identify the 
most relevant collective mode, and its effect on electronic 
excitations in the HTSC. 

In this Letter, we propose that STS experiments can 
directly address this issue by studying the effects in the 
LDOS of a dj.2_ y 2-wave superconductor that arise from 
a coupling between impurities and collective spin and 
charge modes. In particular, we show that such a cou- 
pling induces a static spin (charge) density droplet which 
represents an extended scattering potential for the elec- 
tronic degrees of freedom and thereby affects the su- 
perconductor's local electronic structure. Indeed, static 
spin droplets around Ni and Zn impurities have been 
observed in nuclear magnetic resonance (NMR) experi- 
ments EH, El- By using the T-matrix dEET 
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and 

Bogoliubov-de Gennes (BdG) [l5( formalisms, we demon- 
strate that scattering off spin and charge droplets gives 
rise to qualitatively different fingerprints in the supercon- 
ductor's LDOS. In particular, we find that the coupling 
of an impurity to a spin mode peaked at Q = (tt, it) pre- 
vents the creation of a resonant impurity state, leads to 
characteristic spatial patterns of the spin-resolved LDOS 
directly reflecting the mode's antiferromagnetic nature, 
and gives rise to a pronounced magnetic field dependence 



of the LDOS. Since all of these features are absent in the 
presence of a charge droplet, our results demonstrate that 
impurities can be used to identify the nature and relative 
strength of collective modes not only in the HTSC, but 
in strongly correlated electron systems in general. 

The coupling of a single impurity with spin Si mp lo- 
cated at site R to collective spin and charge modes, rep- 
resented by the operators s e and n e , respectively, is de- 
scribed by the Hamiltonian [l|| E3| 



(1) 



with J, U > and (ni mp ) = 1. Note that the creation of a 
static spin droplet requires a non-zero spin polarization of 
the impurity, (Sf mp ), which can be induced, for example, 
by applying a magnetic field H <C in the a6-plane 
thus avoiding complications arising from the creation of 
vortices 0. Experimentally, (S? mp ) = CH/(T + 6) 
obeys the Curie- Weiss law 12| . The impurity- mode cou- 
pling of Eq.( 
oscillations 



TJ induces static spin [1| 
171 ] described by 



and charge density 



<s*(r)> = J(S! mp )x s (r-K,LO = 0) , 
(Sn e (r)) = U(n imp )xc{r = 0) , 



(2) 



respectively. Here, Xs(Xc) is the spin (charge) suscep- 
tibility, Sn e (r) = n e (r) — no, and no is the uniform 
charge density. For the static susceptibilities in mo- 
mentum space, we make the ansatz \ s . c (q,uj = 0) = 
Xo V(£^c + (ci — Qs,c) 2 ) where f s , e is the respective cor- 
relation length, Q s = (tt, tt) [201 ] and hence = 0, in 
agreement with NMR experiments @, E3, El [12, E3, Ei | . 
The mean-field Hamiltonian of the entire system is given 
by 
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9s(sz{r))cr z a/} - g c (5n(r))l a p] c^c^ , (3) 



where c\ a creates an electron with spin a at site r, t rr i 
is the hopping integral between sites r and r', and A rjI ./ 
is the c? a; 2„ 1) 2-wave superconducting (SC) gap. The last 
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term in Eq.([3]) describes the scattering of electrons by 
the total spin density (s z (r)) = {s z e (r)) + (S? mp )S r ,R 
and effective charge density (Sn(r)) = (Sn e (r)) + 
&vM.{ n imp) ■ The droplets' scattering strength is de- 
termined by only two parameters: for a spin droplet 
by r] s = Jxs(0,0) and g s = g a (S? mp ), such that 
9s(s z {r)) = g s [VsXs(r - R,0)/x s (0,0) +<J R>r ], and for a 
charge droplet by rj c = Uxc(0, 0) and g c = g c (ni mp ). 

We study the effects of electronic scattering on the 
LDOS by using two complementary methods: the T- 
matrix (l3l . approach, which allows us to investigate 
large host systems but assumes a spatially constant SC 
order parameter (SCOP), and the Bogoliubov-de Gcnncs 
(BdG) [f 5] formalism, which accounts for spatial varia- 
tions of the SCOP, but can only treat small system sizes. 
Within the T-matrix approach, the Green's function ma- 
trix in Matsubara space is given by [lj] 

G(r,r',w n ) = G (r,r',w„) 

+ ^Go(r,l,^)r(l,p, W „)Go(p,r', W „) , (4) 
i.p 

where the sum runs over all droplet sites. The T-matrix 
is determined from 



f(l,p,w B ) = Vitfi,p 



fi£6o(l J s ) w n )f(s,p,w n ) , (5) 



when only one of the modes couples to the impurity and 
present in Fig. [lja) [(b)] the normalized spin (charge) 
density for a pure spin (charge) droplet [the center of 
the droplet is located at (0, 0)] [23|. To directly compare 
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where Go = [iu> n ao — £kf3 + At^ip is the unperturbed 
Green's function, V\ = —g s (s z (l))<r z + g c (Sn(l))l, — 
—2t (cos k x + cos k y ) — At' cos k x cos k y — [i is the normal 
state tight binding dispersion with t = 300 meV, t'/t = 
-0.4, and fi/t = -1.083, representative of the HTSC [13], 
and Afc = A (cosfc x — cosfc y )/2 is the SC gap with A = 
30 meV. The LDOS, N(r,w) = A n (r,cj) - A 22 (t,-cj) 
with An (t,u>) = — Im Gu (r,u) + iS)/TT and 6 — 0.2 meV is 
obtained from Eq.(j4]). In contrast, in the BdG formalism 
[lit , one solves the eigenvalue equation 

<•> 

with Hf r , = t rr < + (Tg s (s z (r)) - g c (n(r)) - ^)5 T ^, and 
self-consistently computes the SC gap via 

A rr ' = -— ^ MrK(r') + u n (r')v n (r)] tanh (^j; 

(7) 

where the sum runs over all eigenstates of the system, and 
V = 0.7375 t yields the same (clean) SC gap as taken in 
the T-matrix approach. The LDOS is obtained via 

N(cj, r)=J2 [«n(*)*(" - E n) + v 2 n {v)5(uj + E n )] . 

n 

In order to identify the qualitative effects of spin and 
charge droplets on the LDOS, we first consider the case 



FIG. 1: (color online) Normalized spin (a) and charge (b) 
density along r = (x,0). LDOS inside a spin (c) and charge 
(d) droplet with ^ SlC = 1 and g a ,§c = 400 meV (the dashed 
lines represent the clean LDOS). Evolution of the LDOS with 
increasing rj 3tC for a spin (e) and charge (f) droplet. 

the effects of the droplets, we use g s = g c and £ c = £ s 
with £ s = 5a representative of the underdoped HTSC 
[IH and present in Figs. Q] (c) and (d) the LDOS ob- 
tained from the T-matrix approach inside the spin and 
charge droplet, respectively. The low-frequency LDOS in 
these two droplets exhibits significant qualitative differ- 
ences. Inside the spin droplet, the SC coherence peaks 
are clearly visible and no impurity resonance exists in- 
side the SC gap (the LDOS exhibits, however, weak 
Friedel-like oscillations). In contrast, inside the charge 
droplet, the SC coherence peaks are strongly suppressed 
and a resonant impurity state is formed with peaks in 
the LDOS at ±2 meV. These qualitatively different ef- 
fects of the charge and spin droplet on the LDOS arise 
from the spatial forms of their scattering potentials. For 
a spin droplet, the alternating sign of (s 2 (r)} and hence 
of the scattering potential leads to destructive interfer- 
ence of scattered electrons which prevents the creation 
of an impurity state inside the SC gap. This interpre- 
tation is supported by the fact that the impurity reso- 
nance of the decoupled (r) s = 0) impurity [see Fig. H{e)] 
shifts to higher energies with increasing r) s , implying a 
decrease in the effective scattering strength of the im- 



3 



purity. Note that for sufficiently large rj s , the LDOS is 
qualitatively different from that of a decoupled impurity 
(r) s = 0). In contrast, for the charge droplet, the scat- 
tered electrons interfere constructively due to the same 
sign of the scattering potential at all sites, which leads to 
an increase in the effective scattering strength. Accord- 
ingly, the impurity resonance of the decoupled (rj c = 0) 
impurity shifts to lower energies with increasing rj c [see 
Fig-IHf)]- The modes' different momentum dependence 
thus leads to distinct fingerprints of a spin [Fig. [He)] and 
charge [Fig.QJd)] droplet in the LDOS. 

Another important fingerprint of the spin droplet can 
be found in the spatial structure of the spin-resolved 
LDOS, as shown in Fig.^a) and (b) where we present the 
spin-t and spin-j LDOS at adjacent sites. As expected 
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FIG. 2: (color online) Spin resolved LDOS for g a = 400 meV 
at adjacent sites (a) r = (1,0), and (b) r = (2,0). Intensity 
plot of the (c) spin-j and (d) spin- J. LDOS at lu cp = 30 meV. 



from the spatially alternating sign of the scattering po- 
tential, the frequency dependence of the spin-t LDOS at 
r = (1,0) is qualitatively similar to that of the spin-j 
LDOS at r = (2, 0), and vice versa. This behavior is ob- 
served for all sites of the droplet that belong to different 
(antifcrromagnetic) sublattices. As a result, the spin-j 
and spin-j. LDOS possess a spatially complementary in- 
tensity pattern, as shown in Figs. [He) and (d) where we 
present an intensity plot of the spin-resolved LDOS at 
the frequency of the hole-like coherence peak, uj cp = 30 
meV. At sites where the spin-j LDOS is large, the spin-j 
LDOS is small, and vice versa. In contrast, in a charge 
droplet, the spin-j and spin-j LDOS are identical. Hence, 
the qualitatively different behavior of the spin-resolved 
LDOS in spin and charge droplets is directly linked to 
the magnetic/non- magnetic nature of the modes. 

Another qualitative difference between spin and charge 
droplets arises from the magnetic field dependence of the 



impurity's spin polarization (Sf mp (H,T)) = CH/(T + &) 
[12| , Since g s = g s (S* mp (H,T)), it immediately follows 
that the droplet 's scattering strength and hence the re- 
sulting LDOS are expected to change with H and T. 
To demonstrate this effect, we present in Fig. [SJa) the 
total LDOS for several values of g s representing differ- 
ent magnetic fields. Consider, for example, that a given 
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FIG. 3: (color online) (a) Total LDOS for a spin droplet at 
r = (0,0) for several g s = g s {S! mp (H,T)) (U = 0). (b) Total 
LDOS for a droplet with (s*(r)>, (8n(r)) ^ and a = 1/3. 

Hq <C corresponds to g s = 100 meV. Increasing the 
magnetic field to 2H , (g s = 200 meV) or 4H (g s = 400 
meV) leads to a suppression of the LDOS in the droplet 
[24j ]. Since the coupling between impurity and a charge 
mode is unaffected by the magnetic field, the observation 
of a magnetic field dependent LDOS as shown in Fig.[3ja) 
is another direct signature of a static spin droplet. 

If both charge and a spin modes are present in the su- 
perconductor and simultaneously couple to an impurity, 
we find that properties of the resulting LDOS are pre- 
dominantly determined by the ratio a = g c 7y c /(g s ?y s ), as 
shown in Fig. [21b) for a — 1/3. As a changes from 
to 00, the LDOS changes continuously from being "spin- 
like" [FigOJc)] to being "charge-like" [FigUJd)]. We ex- 
pect that the relative strength of the two modes can be 
determined by studying the iJ-dependence of the LDOS 
and from the form of the spin-resolved LDOS. 

The pairbreaking nature of impurities in d x 2_ y 2--w&ve 
superconductors leads to the suppression of the SCOP, 
an effect which is not taken into account in the T-matrix 
approach. To study the relevance of this suppression, we 
compare the results of the T-matrix approach with those 
of the BdG formalism. In Fig. 0Ja) we present the spa- 
tial form of the SCOP [see Eq.([7])] in a spin droplet for 
a system with N = 60 x 60 sites and the same parame- 
ters as above. The SCOP is significantly suppressed only 
near the center of the droplet, and quickly recovers its 
bulk value within a few lattice spacings from the center. 
In order to ascertain how the SCOP's suppression affects 
the LDOS, we compare in Fig. |4jb) [Fig. Q|c)] the re- 
sults of the T-matrix and BdG approaches for the spin-j 
(spin-j) LDOS at oj cp . The LDOS obtained from both 
approaches is in good qualitative agreement and exhibits 
the same oscillations characteristic of the antiferromag- 
netic nature of the droplet. The quantitative differences 
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FIG. 4: (color online) (a) SCOP along r = (x, 0) for a spin 
droplet. Inset: intensity plot of the SCOP, (b) Spin-| and (c) 
spin- J. LDOS of the T-matrix (red line) and BdG (black line) 
approaches at u) cp = 30 meV along r = (x,0). 



are minimal: first, the LDOS of the BdG approach un- 
dergoes a 7r-phase shift at Ar m 10ao from the center 
of the droplet, which is absent for the T-matrix results. 
Second, the spin-| LDOS of both approaches is out-of- 
phase at r = (0,0). A more detailed analysis shows that 
the latter effect arises from a decrease of the droplet 's ef- 
fective scattering strength due to the suppression of the 
SCOP. Since similar good agreement is obtained for a 
charge droplet, we conclude that the suppression of the 
SCOP has only minor quantitative effects on the LDOS. 

While we considered above a static spin droplet, recent 
studies [25l | suggest that our results remain valid even for 
H = as long as the fluctuation time of the droplet sat- 
isfies ts S> 1/Ep. If the spin excitations are sufficie ntly 
damped, ts — > 00 and the droplet becomes static [261 ] - 
One might therefore wonder whether signatures of an im- 
purity coupling to a quasi-static mode have already been 
observed [27|. The effects of collective modes on quasi- 
particle interference patterns in the fourier-transformed 
LDOS were discussed in Refs. 
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Finally, placing an impurity into the Cu02-plane of the 
HTSC can lead to static lattice distortions. Since such 
a "phonon droplet" is non-magnetic in nature, it gives 
rise to a scattering term in Eq.fl]) that is identical to 
that of the charge droplet albeit with a potentially more 
complicated momentum dependence of g. We therefore 
expect the qualitative effects of static phonon and charge 
droplets to be quite similar. 

In conclusion, we have shown that impurity induced 
static spin or charge droplets exert qualitatively differ- 
ent effects on the LDOS of a d x 2_ y 2 superconductor. We 
find that in a spin droplet, the creation of a resonant 
impurity state is suppressed, the spin-resolved LDOS ex- 
hibits complementary spatial patterns, and a character- 
istic magnetic field dependence. Since these fingerprints 
are absent in a charge droplet, our results demonstrate 
that impurities can be used to identify the nature and rel- 
ative strength of collective modes not only in the HTSC, 
but in strongly correlated electron systems in general. 
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